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We consider a multicomponent atomic Bose-Einstein condensate optically trapped in a far-oflt 
resonant dipole trap. Drawing an analogy with the optical situation, we show that this system can 
be regarded as an matter-wave analog of optical multiwave mixing. We concentrate specifically on 
condenstates in the hyperfine ground state _F = f , in which case a simple analogy with optical four- 
wave mixing can be established. This opens up the way to realize matter-wave phase conjugation, 
whereby an atomic beam can be "time-reversed." In addition to transfering population between a 
"central" mode and incident and retroreflecting beams, matter- wave phase conjugation also offers 
novel diagnostic tools to study the coherence properties of condensates, as well as to measure the 
relative scattering lengths of hyperfine sublevels. 

PACS numbers: 03.75.-b 03.75.Fi 05.30.Jp 42.65.Hw 
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I. INTRODUCTION 

The experimental observation of Bose-Einstein con- 
densation in low density atomic vapors |l| has triggered 
a flurry of theoretical activity Theoretical predic- 
tions on the condensate dynamics, ground state popula- 
tion and spectrum of elementary excitations have been 
made and are in excellent agreement with experiments. 

At the same time, further experimental advances have 
led to the realization of multicomponent condensates, 
both in ^^Rb and in ^'^Na. In the first case, sympa- 
thetic cooling, together with a fortuitous coincidence in 
the scattering lengths of the spin states m = —1/2 and 
m = —3/2 led to the coexistence of both components 
in a magnetic trap Q. A multicomponent condensate 
was also achieved with the three hyperfine ground state 
components of sodium in a far-off resonant dipole trap 
Q. These results, along with further experiments in- 
volving condensates in double well potentials have 
now led to considerable theoretical work on the static 
and dynamic properties of multicomponent condensates, 
including studies of their true ground state analy- 
ses of the elementary excitations spectrum, and the de- 
termination of their instability regions |^]. In addition, 
multicomponent condensates open up the way to novel 
schemes to launch vortices and permanent currents ||] 
and to the study of novel phenomena such as vector and 
quadrupole spin wave modes, the topological and ener- 
getical instabilities of doubly quantized singular vortices 
1^, etc. The two-body interactions characteristic of spin- 
1 condensates can lead to intercomponent coupling via 
two-body collisions. 

In the zero-temperature limit, a g-component Bose 
Einstein condensate can be thought of as a g-mode sys- 
tem, whereby the various modes are coupled by two-body 
(and possibly higher-order) collisions which result in the 
exchange of particles between these modes. As such, they 
correspond to a situation quite similar to that of multi- 
wave mixing in nonlinear optics. A main difference is of 



course that in the case of matter waves, the coupling is 
due to the collisions, which find their origin in the electro- 
magnetic interaction between the atoms. Collisions can 
then be thought of as the effective atom-atom interaction 
resulting from the elimination of the electromagnetic field 
from the system dynamics. This is to be contrasted with 
the optical case, where multimode mixing relies on the 
interaction of the electromagnetic field with a common 
atomic sample whose dynamics is traced over. In that 
case, it is the elimination of the material dynamics that 
results in an effective field- field coupling. This observa- 
tion is the origin of nonlinear atom optics, which is the 
matter- wave equivalent of nonlinear optics. 

A close analogy can easily be established between the 
dynamics of a spin-1 condensate as realized in sodium 
experiments and the situation of degenerate four-wave 
mixing in optics, as we demonstrate explicitly in this pa- 
per. In particular, for situations where the to = state 
is macroscopically populated while the m = ±1 states 
are weakly excited, one can think of the first state as 
a "pump" or "central" mode, while m = ±1 form side 
modes, which are coupled via the pump, leading to the 
familiar effects of degenerate four- wave mixing, including 
phase conjugation. This is the effect that we study in de- 
tail in this work. We then show how matter-wave phase 
conjugation can be used as a diagnostic tool to study the 
coherence properties of Schrodinger fields, as well as the 
relative scattering lengths of the states involved. 

Matter-wave phase conjugation has previously been 
studied, but in a situation where the coupling between 
the partial matter waves was induced by the near- 
resonant electric dipole-dipole interaction As such, 
it relied explicitly on having a substantial population of 
electronically excited atoms, and the incoherent effects 
of spontaneous emission rapidly destroyed the coherent 
wave coupling responsible for phase conjugation. In con- 
trast, the situation with a condensate of dipole-trapped 
sodium atoms does not suffer from this drawback: since 
we are considering ground-state atoms in a far-off res- 
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onant trap with hyperfine levels coupled primarily via 
ground state collisions, spontaneous emission is certainly 
negligible. In addition, the fact that the atoms are in 
a trap changes the situation somewhat from the free- 
space geometry considered in our earlier work, since the 
atomic sample can easily be tightly confined in the trans- 
verse dimensions and hence does not suffer from free- 
space diffraction. Our main result is to demonstrate 
that a trapped condensate can then be used as a phase- 
conjugate mirror for a weak atomic beam, thereby effec- 
tively "time-reversing" it. 

Section II describes our physical model, and derives 
the coupled-wave equations for the three components of 
the condensate in the Hartree regime. This is applied in 
section III to the discussion of matter- wave phase conju- 
gation in two-dimensional atomic traps. We concentrate 
explicitly on the undepleted pump regime, and show how 
the phase conjugate signal depends explicitly on the rel- 
ative scattering lengths of the hyperfine levels involved. 
Finally, the possible experimental verification of our pre- 
dictions, as well as a summary and outlook are given in 
section IV. 



II. PHYSICAL MODEL 

We consider a condensate of ^"^Na atoms in their F = 1 
hyperfine ground state, with three internal atomic states 
|F = l,m = -1), |F = 1,TO = 0} and |F = l,m = 1} of 
degenerate energies in the absence of magnetic fields. It 
is described by the three-component vector Schrodinger 
field 

*(r,i) = {vI/_i(r,t),vI/o(r,i),vI.i(r,i)} (1) 

which satisfies the bosonic commutation relations 

[Mr,t),^l{r',t)] = 5,,S{r-r'). (2) 

Accounting for the possibility of two-body collisions, its 
dynamics is described by the second-quantized Hamilto- 
nian 



n = J dr^'<{r,t)Ha^{r,t) 

+ J{dr}^^{ri,t)<if\r2,t)V{ri-r2)^{r2,t)^{ri,t), (3) 
where the single-particle Hamiltonian is 

Ho = pV2M + Vtrap (4) 

and the trap potential is of the general form 
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Vtra 



U{r)\F = l,m){F = 1, 



(5) 



m— — 1 



Here p is the center-of-mass momentum of the atoms of 
mass M and U{y), the effective dipole trap potential for 



atoms in the |l,m) hyperfine state, is independent of m 
for a non-magnetic trap. 

The general form of the two-body interaction V{yi — 
has been discussed in detail in Refs. |p|,pl|. We repro- 
duce its main features for the sake of clarity. Consider 
situations where the hyperfine spin i^j = 1 of the indi- 
vidual atoms is preserved. We label the hyperfine states 
of the combined system of hyperfine spin F = Fi + F2 
by \f,m) with / = 0,1,2 and m = —/,...,/. In the 
shapeless approximation, it can then be shown that the 
two-body interaction is of the general form H 



where 



V{ri - V2) = (5(ri - ra) ^ Tiff/T'/, 



gf = Anhaf /M, 



(6) 



(7) 



Vf = ^„ I/, m)(/, m| is the projection operator which 
projects the pair of atoms into a total hyperfine / state 
and af is the s-wave scattering length for the channel 
of total hyperfine spin /. For bosonic atoms only even 
/—states contribute, so that 



V{ri - r2) = M(ri - r2){g2'P2 + goVo) 

= ■^'^(ri - ^2) (co + C2F1 • F2) 



(8) 



In this expression, 



Co = 2(.go + 2.92)73 
C2 = 2(^2 - .go)/3. 



(9) 



Substituting this form of V{ri — into the second- 
quantized Hamiltonian (b|) leads to 



P 

2M 



C/(r) 



*m(r,t) 



7i f 

+ - / dr{(co +C2)[«'I^'I*i*i + 

-f2(co - C2)«'l*l*Li*-l 

-|-2c2(*I*^i*o*o + ff.c.)}. (10) 

This form of the Hamiltonian is quite familiar in quan- 
tum optics, where it describes four- wave mixing between 
a pump beam and two side- modes, which are identified 
with the field operators ^0 and ^'±i in the present sit- 
uation. Specifically, we observe that the three terms in 
the two-body Hamiltonian which are quartic in one of 
the field operators only, i.e. of the form ^'^^^I'l^I'i^I'i can 
be readily interpreted as self-defocussing terms, corre- 
sponding to the fact that the two-body potential is, for 
a positive scattering length and a scalar field, analogous 
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to a defocussing cubic nonlinearity in optics. The terms 
involving two "modes", i.e. of the type ^'■VE'iVE'jvE'j, con- 
serve the individual mode populations of the modes and 
simply lead to phase shifts. Finally, the terms involving 
the central mode ^'o and both side-modes are the con- 
tributions of interest to us, since they correspond to a 
redistribution of atoms between the "pump" mode \l/o 
and the side-modes V^ii, e.g. by annihilating two atoms 
in the central mode and creating one atom each in the 
side-modes. This is the kind of interaction that leads to 
phase conjugation in quantum optics, except that in that 
case the modes in question are modes of the Maxwell field 
instead of the Schrodinger field. Note also that a similar 
mechanism is at the origin of amplification in the Collec- 
tive Atom Recoil Laser (CARL) |l2| , |l3| , except that in 
that latter case, the Schrodinger field mode coupling is 
induced by optical transitions in the atoms. 

In the Hartree approximation, which is well justi- 
fied for condensates at T = 0, the many-body prob- 
lem reduces to an effective single particle problem for 
the Hartree wave function (/)„j(r,t). It is easily shown 
that to be governed by the system of coupled nonlinear 
Schrodinger equations 



i0_i(r,i) = i 



— + U(r) 
2M ^ ' 



+ [(co + C2)(|0-ir + 100 I ) + (co - C2) 
^2 



i(po{r,t) 



2M 



U{r) 



TV {C2 00 00 01 

0o + A^{co|0ol'0o 



+ (C0 + C2)(|0-lP + \(f>l\^)(t>0 + 2C2010_10^} 



i0i(r,t) = - 



01 + 7V{C20O0O0*1 



+ [(co + C2)(|01p + |0OP) + (co - C2)| 



l}. (11) 



Just as in the familiar quantum optics case, we consider 
in the following a situation where the central mode, de- 
scribed by the Hartree wave function 0o , is strongly pop- 
ulated initially, while the side-modes 0±i are weakly pop- 
ulated. In other words, we consider the phase conjuga- 
tion of a weak atomic beam from a reasonably large con- 
densate. In that case, it is appropriate to introduce the 
matter-wave optics equivalent of the undepleted pump 
approximation, whereby 



0. 



(12) 



In that case, the problem reduces to a set of coupled 
mode equations for the two side-modes 0±i, the central 
mode acting as a catalyst for the coupling between them. 



III. PHASE CONJUGATION IN DIPOLE TRAPS 

In what follows we consider atomic samples confined 
in two-dimensional harmonic trap. The trap potential 



U{r) which is as we recall independent of the atomic in- 
ternal state m Q for a dipole trap, is taken to be of the 
harmonic form 



U{r) 



y^)/2 



(13) 



for simplicity. That is, we assume that the dipole trap 
confines the atoms in the transverse plane {x,y), but 
not in the longitudinal direction z. This geometry al- 
lows one to consider side-modes propagating along that 
axis, rather than bouncing back and forth in an elongated 
trap. In case of tight confinement in the transverse di- 
rection, we can assume to a good approximation that the 
transverse structure of the condensate is not significantly 
altered by many-body interactions and is determined as 
the ground-state solution of the transverse potential. 

Expressing the Hartree wave function associated with 
the hyperfine level m as 



(j)mir,t) ^ ip^{x,y)ipm{z,t)e 



(14) 



we then have 
huJoip±{x,y) 



\-- 
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2M 




^ dy^J 



Muj?, 



(p±{x,y). (15) 



Substituting this expression into Eqs. (|11|) and project- 
ing out the transverse part of the wave function yields 
the coupled one-dimensional Gross-Pitaevskii (coupled- 
mode) equations 

+ [(C0 +C2)(|v3_i|2 + |(^on + (co - C2)|(pi|2]^_i} 
+ (co + C2)(|(^_iP + \ipi\^)ipo + 2c2ipi^P-lLpl} 

iipi{z,t) = ~^-Q^(pi + Nr]{c2(po(pQip*_i 

+ [(C0 + C2)(|^lP + \ipo\^) + (Co - C2)|^_lP]^i}, (16) 



where 



V 



J dxdy\(p±{x,y)\'^ 
J dxdy\(p±{x,y)\'' 



(17) 



The physical situation we have in mind is that of a weak 
"probe" in the hyperfine state m = — 1 propagating to- 
ward a large condensate in state m = and at rest in the 
dipole trap, and generating a backward-propagating con- 
jugate matter wave in the hyperfine state m = — 1 (see 
Fig. 1). Hence we express the longitudinal component of 
the Hartree wave function as 



<p(z,t) 



ip-iiz,t) 
ipo{z,t) 
(fi{z,t) 



^_i(z,t)e-^*=^e-"' 
2V'o cos(A:z)e-*'^* 
^pl{z,t)e 



ikz ^—iut 
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FIG. 1. Geometry of the matter wave phase conjugation 
with 3-component vector fields. 

where the slowly varying envelopes of the Hartree 
wave function components to = ±1 satisfy the familiar 
inequalities 



52 Q 
I— ^Vm| < k\—i>J\ < fc^lV'n 

az^ oz 



(18) 



and we have additionally invoked the undepleted pump 
approximation (12). Note that in this ansatz the "pump" 
wave function (^0 is described by a standing wave. This 
spatial structure is required in order to achieve momen- 
tum conservation, a direct consequence of the fact that 
a standing wave can be viewed as a superposition of two 
counterpropagating atomic waves. A state with such a 
periodic spatial structure can be achieved for instance by 
interfering two condensates ||l5| , in a grating matter- wave 
interferometer, jl^ or in CARL ||l^. To the first order 
in the probe and signal fields, this geometry leads to a 
linearized system of two coupled-mode equations for the 
probe and condensate fields. In the stationary state they 
reduce to 

fik d 



2M dz 



hk d 

'2Md~z'^*^^'^ " -^^'?[2(co + C2)poV't(^) 



C2i^li^liz)], 

+ C2lplij-l{z)], 



(19) 



where po = l^oP- 

The form of these equations is familiar from optical 
phase conjugation and their solution is well-known. Be- 
fore giving them explicitly, though, we note that they 
contain two contributions. For instance, the equation 
for the phase conjugate wave ip* contains a term pro- 
portional to the density po of the condensate and the 
field itself. In the absence of the second term, it would 
simply lead to a phase shift of ipi- Physically, it results 
from the self-interaction of the conjugate field, catalyzed 
by the condensate (pump) component. Its origin can be 
traced back to the term proportional to ^J^'J^'o^'o in the 
Hamiltonian ([l0|). The second term, in contrast, couples 
the two side-modes via the condensate and is responsible 
for phase conjugation. Note that it is not proportional to 
the condensate density poj but rather to ipQ. We return 
to this point later on. 

The general solution of Eqs. ( |l9|) reads ||l^ 



V'-i(z) 



( - te""^ sm{\K\z)i;*{L) 



cos(|k|L) 
+ cos(|k|(z-L))V'_i(0)) 

( cos(|K|z)i/'i(i) 



cos(|k|L) 

+ i e-'^^ sm{\K\{z ~ L))^j*_^{0)) , (20) 



where 



a = 2Nri{co + C2)po, 



Nr]C2ipo 



nk/2M ' 



and 



J0 



k/\k\. 



(21) 



(22) 



(23) 



For the probe ■(/'_i(0) incident at z = and no incom- 
ing conjugate signal ipiiL) = 0, the conjugate wave in 
the input plane z = becomes 



^i(O) = -ie-"3tan(|K|L)V'*i(0), 



(24) 



which demonstrates that the interaction of the probe and 
the condensate results in the generation of a counterprop- 
agating phase-conjugated signal. Note that the intensity 
of the conjugate wave exceeds that of the incoming wave 
for 7r/4 < \k\L < 3tt/A and phase conjugation oscilla- 
tions (PCO) ||l^ can occur for \k\L = tt/2, the so-called 
oscillation condition. 



IV. EXPERIMENTAL FEASIBILITY AND 
OUTLOOK 

In order to determine the feasibility of matter-wave 
phase conjugation in state of the art experiments, we 
briefly discuss the values of the oscillation parameter \k\L 
that can be achieved in current ^"^Na BEC experiments. 



From the definition (22) we have 



\k\L - Ni] 



Q2 - Qq 
3fc ' 



(25) 



where we have taken that due to normalization "ipgL 
PqL ^ 1 and that |^ 



C2 = 47rS,(a2 - ao)/3M 



(26) 



with oo and 02 being the singlet and triplet state scat- 
tering lengths respectively. For sodium, these scatter- 
ing lengths are estimated as ^ (02 — ao)/3 '~ 0.04a2 
lO^^^m. In the MIT optical confinement experiments 
the number of trapped atoms is of the order 5 — 10-10 
and the transverse dipole trap frequency loq is of the order 
of 10* sec~^, so that the transverse ground state size of 
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0.25/im. FromEq.ll 



As a result the oscillation 



the condensate aj 
we have 77 ^ aj^ ^ 10^'^r 
parameter is \k\L ^ 10^"/fc where fc is as we recall the 
wave number of a Schrodinger field. In case the conden- 
sate sidemodes are obtained by diffraction on a standing 
light wave jl6j, we have k — 27r/A ~ lO^m^^ and thus 
\k\L ^ 10'^. This means that the oscillation condition 
\k\L = 7r(2n + l)/2 (n - integer) can be met in current 
BEC experiments. 

In addition to its interest from a nonlinear atom optics 
point of view, matter-wave phase conjugation could also 
be used as a diagnostic tool for Bose-Einstein conden- 
sates. For instance, we noted that the parameter \k\L 
is proportional to the difference in scattering lengths be- 
tween the singlet and triplet states. Hence, this quan- 
tity could in principle be inferred from phase conjugation 
measurements. In addition, we recall that the phase con- 
jugate signal is not determined by the condensate density 
Pq, but rather by ■0o- While the distinction between the 
two is expected to be minimal for large condensates, and 
is essentially ignored in the Hartree and undepleted pump 
approach of the present paper, this will no longer be the 
case for smaller condensates. In such situations, phase 
conjugation provides one with a probe of the coherence 
properties of the condensate. Future work will analyze 
these aspects of the problem, as well as the role of higher- 
order correlation functions in the atom statistics of the 
phase-conjugate mode. 
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